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Abstract: Black holes in asymptotically Lifshitz spacetime provide a window onto 
finite temperature effects in strongly coupled Lifshitz models. We add a Maxwell 
t> ■ gauge field and charged matter to a recently proposed gravity dual of 2+1 dimen- 

^ sional Lifshitz theory. This gives rise to charged black holes with scalar hair, which 

correspond to the superconducting phase of holographic superconductors with z > 1 
Lifshitz scaling. Along the way we analyze the global geometry of static, asymptot- 
ically Lifshitz black holes at arbitrary critical exponent z > 1. In all known exact 
solutions there is a null curvature singularity in the black hole region, and, by a 
general argument, the same applies to generic Lifshitz black holes. 
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1. Introduction 

A holographic superconductor [1, 2] is dual to an asymptotically AdS spacetime with 
a charged black hole carrying scalar hair. It is a theoretical construction where a 
superconducting phase transition in a strongly coupled system is studied via the 
AdS/CFT correspondence [3]. Having a black hole places the system at finite tem- 
perature, the black hole charge gives rise to a chemical potential in the dual theory, 
and the scalar field hair signals the condensation of a charged operator in the dual 
theory. Without a chemical potential all temperatures are equivalent due to the 
underlying conformal symmetry and there can be no phase transition. With a chem- 
ical potential, on the other hand, a new scale is introduced that allows for a phase 
transition at some critical temperature. 
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The black holes in question are found as solutions of AdS gravity coupled to a 
Maxwell gauge field and matter in the form of a charged scalar. Black holes that 
are neutral under the Maxwell field do not develop any hair. Near a charged black 
hole, however, the gauge field sourced by the black hole couples to the charged 
scalar and induces a negative mass squared sufficient for condensation provided the 
temperature is low enough. It was observed in [4] that even a neutral scalar can 
lead to condensation below a non zero critical temperature due to the existence of 
a new and effectively lower Breitenlohner-Freedman bound [5] near the horizon of a 
low temperature, near extremal, AdS-Reissner-Nordstrom black hole. See [6] and [7] 
for recent reviews of holographic superconductivity. 

In the present paper we show that similar phenomena can occur in a theory 
exhibiting Lifshitz scaling, 

t^XH, x^Ax, (1.1) 

with z 1. Models with scaling of this type have, for example, been used to 
model quantum critical behavior in strongly correlated electron systems [8]- [11]. Our 
starting point is the model proposed by Kachru, Liu, and Mulligan [12] for the 
holographic study of strongly coupled 2+1 dimensional systems with Lifshitz scaling. 
While this model in itself does not support a phase transition to a superconductor 
it turns out that a relatively simple extension does. 

A quantum critical point exhibiting dynamical scaling of the form (1.1) has a 
gravitational dual description in terms of a spacetime metric of the form 

ds 2 = L 2 (-r 2z dt 2 + ^ + r 2 d 2 ^j , (1.2) 

which is invariant under the transformation 

t^\ z t, r->^, x^Ax. (1.3) 
A 

Here L is a characteristic length scale and the coordinates (t, r, x 1 , x 2 ) are taken to 
be dimensionless. Critical points with z > 1 are often said to be non-relativistic. 
This can be motivated by considering null geodesies at fixed r = r in the Lifshitz 
background (1.2). Along such a geodesic one finds that (dx/dt) 2 = r 22 " 2 , and so for 
z > 1 the effective speed of light in the boundary theory diverges as r — > oo. 

In [12] it was shown how the fixed point geometry (1.2), with z > 1, can be ob- 
tained from an action coupling 3+1 dimensional gravity with a negative cosmological 
constant to abelian two- and three-form field strengths, 

S = J d A x^j (R-2A) (1.4) 
-\\ *F {2) A F (2) - X - J *H {3) A H {3) - c I B {2) A F (2) 
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with Hrg\ = dB(2) and the cosmological constant and the coupling between the form- 
fields given by A = — z2 ^ 2 +4 and c = The equations of motion for the form 
fields can be written 

d * F(2) = -cH(3), d * if (3 ) = -cF( 2 ), (1.5) 
and the Einstein equations are 

Gfiv + Ag^u = \{F»xF v x - - A g, u F Xa F x °) + \{H^ C H X ° - l -g, u H Xap H x ^). (1.6) 

Black holes in this 3+1 dimensional gravity theory were considered in [13], where 
numerical black hole solutions were found at z = 2 and used to study finite temper- 
ature effects in the dual 2+1 dimensional system. Related work on Lifshitz black 
holes can be found in [14], where topological black holes with hyperbolic horizons 
were included, and in [15, 16], where black holes at general values of z were consid- 
ered. 1 These Lifshitz black holes carry a charge that couples to the two- form field 
strength Ft 2 ), but they are in many respects more analogous to AdS-Schwarzschild 
black holes than charged AdS-Reissner-Nordstrom black holes. In particular, since 
the black hole charge cannot be varied independently of the black hole area, there is 
only a one-parameter family of black hole solutions for a given value of z. 

We will see below that without additional ingredients the underlying Lifshitz 
symmetry prevents the system from undergoing phase transitions. We therefore 
extend the system to include a Maxwell field, with field strength J-(2), and a scalar 
field ip that is charged under the new gauge field but neutral under the original 
Lifshitz fields and H^y This turns out to be sufficient in order to observe 
a superconducting phase transition characterized by Lifshitz scaling with z > 1. 
In our modified theory it is the Maxwell field J r ( 2 ) that corresponds to physical 
electromagnetism while the original Lifshitz gauge fields are viewed as an auxiliary 
construction, whose only role is to modify the asymptotic symmetry of the geometry 
from AdS to Lifshitz. 

The plan of the rest of the paper is as follows. We begin in Section 2 with a brief 
review of Lifshitz black holes in the model (1.4). We extend previous treatments 
by considering the global geometry, including the black hole interior. In those cases 
where an exact solution is known, we find a null curvature singularity at r = and 
a Carter- Penrose diagram as shown in Figure 1. We show that a null singularity is 
generic for black holes in this model. In Section 3 we generalize the model to include 
an additional Maxwell field. We exhibit a family of exact solutions at z — 4, which 
describe Lifshitz black holes that are charged under the new Maxwell field and obtain 
numerical solutions for charged black holes at other values of z > 1. These black 

1 Black hole solutions in other gravity models exhibiting Lifshitz scaling have been considered in 
[17, 18, 19]. 
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holes can be viewed as the Lifshitz analog of AdS-Reissner-Nordstrom black holes. 
In Section 4 we further extend the model by adding a charged scalar field and look 
for the black hole instability that signals the onset of superconductivity. Finally, we 
wrap up with some final remarks in Section 5. 



2. Lifshitz black holes revisited 

The action (1.4) is known to have spherically symmetric static black hole solutions 
of the form 

ds 2 = L 2 (^-r 2z f (rf dt 2 + 9 -^fdr 2 + r 2 (d9 2 + X (0) 2 dip 2 ) j , (2.1) 

with 

{sin6 if k — 1, 
e if k = o, 
sinh# if k = -l, 

where k = +1,0,-1 corresponds to a spherical, flat, and hyperbolic horizon re- 
spectively. An asymptotically Lifshitz black hole with a non-degenerate horizon has 
/(r), g(r) — > 1 as r — y oo and a simple zero of both f(r) 2 and g(r)~ 2 at the horizon 
r = r . 

2.1 Exact solutions and global extensions 

Several families of numerical solutions of this general form have been found [13, 14, 
15] along with a couple of isolated exact black hole solutions. These are a z = 2 
topological black hole with a hyperbolic horizon [14], 



w-serV 1 -^' (2 - 2) 



and a z = A black hole [15], 



k 3k 2 

/(r) = W) = V 1 + 10^ " 400^' (2 ' 3) 

With k — ±1. The k — +1 case has a spherical horizon while k — — 1 corresponds 
to a topological black hole with a hyperbolic horizon. When we extend the system 
to include an additional Maxwell gauge field (see Section 3 below) we will see that 
the exact z = 4 solution (2.3) is a special case of a one-parameter family of exact 
solutions in the extended theory. 

We start our discussion by working out the global geometry of the exact z = 2 
topological black hole (2.2) and constructing the associated conformal diagram. In 
order to study the interior geometry of a black hole one looks for coordinates that 
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are non-singular at the black hole horizon and allow the solution to be extended into 
the black hole. The first step is to transform to a tortoise coordinate r* for which 



ds 2 



21 r 1 :i ^) - dr 2 ) ..- r 2 (>/\ 2 4 sink 2 V ^ 2 ) 



(2.4) 



Here r is to be viewed as a function of r* determined via 



dr 
dr* 



r 3 I 1 



1 

2^2 



which integrates to 



1 



r = 



(2.5) 



(2.6) 



v 2yl — exp(r\,, — rj°) 

As usual, the tortoise coordinate r* goes to — oo at the horizon but note that it goes 
to a finite value in the r — > oo asymptotic region. This is easily seen to be a generic 
feature of asymptotically Lifshitz black holes for any z > 1. 

Next we form the null combinations v = t + r*, 
■u = t — r* and perform a conformal reparametrization, 



(2.7) 



v — >■ V — exp 

■u — > [/ = — exp 

which brings the metric into the form 

ds 2 = L 2 {-Ar A dU dV + r 2 {d6 2 + sinh 2 6d V 2 )) . (2.8) 

In these coordinates the geometry is manifestly non- 
singular at the event horizon and the solution can be 
extended to r < 1/ \/2. 

In general, when extending a black hole solution 
through a non-degenerate horizon the required confor- 
mal reparametrization is of the form 




v — > V = a exp (kv), 
u — )> U = —a exp (—/t-u), 



Figure 1: Conformal dia- 
gram for a Lifshitz black hole. 
The thick lines denote null 
(2.9) curvature singularities. 



where a is an arbitrary constant and k is the surface 

gravity of the black hole in question. The surface gravity determines the Hawking 
temperature, T H = ^, and it is easily checked by other means that T H = is the 
correct value for the z = 2 black hole in (2.2). We have chosen the value of the 
constant a in (2.7) such that UV — > — 1 in the r — > oo asymptotic limit. 
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The z = 2 topological black hole has a curvature singularity at r = 0, which can 
for example be seen by computing the Ricci scalar, 2 

it> = J^(l-22r 2 ). (2.10) 

It follows from 

~ UV = 1 -^ ( 2 - n ) 

that r — > oo corresponds to UV — > — 1 and r — > to UV — > oo. The Carter- Penrose 
conformal diagram 3 in Figure 1 is then obtained by writing 

ix P irQ 
V = tan , U = tan—. (2.12) 

The z = 4 case can be dealt with in a similar manner. There is a null curvature 
singularity at r = and the 

The tortoise coordinate for a z = 4 black hole is given by 

- " r ' = log 4) + log 4) • < 2 ' 13 > 

with fei = ^(2 1 k | — k) and 6 2 = ^(2|A;| +k). The change of coordinates to 

V = exp [&!(&! + b 2 )(v - O] , U = - exp [-^(^ + 6 2 )( M + O] , (2.14) 
renders the metric nonsingular at the horizon, 

ds 2 = L 2 ^-^ ^1 + ^ ^ dUdV + r 2 (rf^ 2 + X \0) dip 2 ) j . (2.15) 

The Hawking temperature of these z = 4 black holes is 

T " = ^ 2 - k ^ < 216 » 

and the Ricci scalar is given by 

1 / 9k 2 3k r \ 

There is a null curvature singularity at r — >• and the conformal diagram is that of 
Figure 1. In Section 2.3 below, we present a general argument that this conformal 
diagram applies to all black hole solutions of the system of equations (1.5) - (1.6). 



2 The fixed point geometry (1.2) is also singular at r = but in a relatively mild way. All scalar 
invariants constructed from the Ricmann tensor are finite but tidal effects between neighboring null 
geodesies diverge as r — >• 0. See f.ex. [6]. 

3 We refer to Figure 1 as a conformal diagram but it should be kept in mind that the boundary 
at r — >• oo is in fact not conformally flat due to the asymmetric scaling of the spatial and time 
directions in asymptotically Lifshitz spacetime. 



6 



2.2 Lifshitz black holes at general z > 1 

The exact solutions offer a glimpse at the parameter space of Lifshitz black holes 
at isolated points. More generic solutions can be obtained numerically. We are 
interested in the global geometry including the black hole interior. For this we take 
our cue from the extension of the exact solution described above and write the metric 
in the form 



The notation is descended from two-dimensional gravity and allows the field equa- 
tions to be witten in a relatively economical way. With this ansatz the metric is 
characterized by two field variables, p and 0, which determine the local scale of the 
v, u plane and the scale of the transverse two- manifold respectively. The null coor- 
dinates v, u are related to the t, r coordinates in (2.1) via v = t + r*, u = t — r*, with 
the tortoise coordinate r* given by 



The Lifshitz fixed point geometry (1.2) has f(r) = g(r) = 1 and in that case 



zr z 

More generally, we see from (2.19) that r\ goes to a finite value asymptotically for 
any Lifshitz spacetime, for which f(r),g(r) — >■ 1 as r — >■ oo. 

2.2.1 The equations of motion 

With the following ansatz for the form fields, 



rfs 2 = L 2 [- e Mv,u) dvdu + e -2*(«,«) (^2 + x (e) 2 d<f 2 )] . 



(2.18) 




(2.19) 
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(2.20) 



F( 2 ) = L f vu (v, u) dv A du, 

i7( 3 ) = L 2 (h v (v, u)dv — h u (v, u)du) A d6 A x(@)d(p, 



(2.21) 
(2.22) 



their field equations (1.5) become 




(2.23) 



(2.24) 



d v {e 2 *h u ) + d u (e 2< X) 




(2.25) 



These can in turn be re-expressed as a single second-order equation, 



d v dj + d v <f>dj + d u <j)d v f + -e 2p f = 0, 



(2.26) 
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where we have defined 

/ = e- 2p ~^f vu . (2.27) 
The equations of motion for and p are obtained from (1.6), 

-d v d u( f> + 2d v <j>d u cf> + -/"^ + (z2+ g +4) e 2 ^ - = 0, (2.28) 

-d v d u p + d v <pd u <P + -e 2 ^ - Ipe 2 " 4 * + ^djdj = 0. (2.29) 
4 Z Zz 

The conformal reparametrization (2.9) will render the metric non-degerate at the 
horizon. The form of the field equations remains the same in the V, U coordinate 
system with a transformed conformal factor, 



2p(«,«) = (- K 2 UV) e 2p{v > u \ (2.30) 



e 

We now adapt a simple numerical method, which was originally developed for the 
study of two-dimensional black holes [20], to the case at hand. We introduce a new 
spatial variable, 

s ee -UV, (2.31) 
and restrict our attention to static configurations. The field equations become 

= sf" + /' + 2s<j>'? - l~fe 2p , (2.32) 
= S <P" + 0' - 2s0' 2 + ^e 2 ^ + i£!±f±!! e 2p _ ife 2 ^, (2.33) 
= sp" + (/- s<P' 2 + -e 2 ^ 2 * - \] 2 e 2 ^ - ^e 4 ^/' 2 , (2.34) 

where primes denote derivatives with respect to s. 
The Lifshitz geometry (1.2) is given by 



/ = y ^-ilr 2 , e"* = r, K 2 s e 2 ^ = r 2 *, (2.35) 

with r(s) obtained from (2.20) and (2.31). The horizon at r = is singular and as 
a result k can take any value in the Lifshitz geometry. The z = 2 topological black 
hole has / = r 2 with r(s) = 1/ a/2(1 — s), while the z = 4 black holes turn out to 
have / = VQ(r 2 + ±), with r(s) obtained via £ = r 5 (l + ^ - ^) and (2.31). 

2.2.2 Numerical solutions 

The event horizon is at s = and s is negative inside the black hole. With the 
assumption of a regular horizon we can read off the following relations among initial 
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values at s = 0, 



0'(O) = I (^-A;e 2 ^ - (Z + ^ +4) + /(0) 2 e^°)J e 2 *°\ (2.36) 
/>) = f/(0)e 2 ' (0) , (2.37) 
p'( ) = 1 ^-^e 2 ^°) + /(0) 2 e 4< ^ e 2 '(°>. (2.38) 

These initial values can now be used to start a numerical integration of the system 
(2.32)-(2.34) from near s = 0, either outwards towards s > or into the black hole 
interior at s < 0. Inequivalent solutions are parametrized by 0(0) and /(0). The 
0(0) initial value gives the value of the area coordinate at the horizon through the 
relation r = e~^, and /(0) determines the magnitude of the radial two-form field at 
the horizon. A shift of p(0) amounts to a global rescaling of the s coordinate and 
does not affect the geometry. 

As discussed in [13], 0(0) and /(0) cannot be varied independently while pre- 
serving the asymptotically Lifshitz character of the geometry. In the z = 2 case 
considered in that paper, there is a zero mode of the linearized system of equations 
near the Lifshitz fixed point, which must be set to zero for the system to approach 
the fixed point geometry (1.2) as r — > oo. This requires a fine-tuning of parameters 
which uniquely determines /(0) in terms of 0(0), or vice versa. In [15] it was shown 
that at z > 2 the zero mode is replaced by a power-law growing mode which must 
be set to zero to obtain an asymptotically Lifshitz geometry. The question was more 
subtle at z < 2 where instead there is a mode with a weak power-law fall-off but a 
finite energy argument [15, 21] may be employed to conclude that this mode must 
also be set to zero. 4 

With the variables that we are using here there is a simple criterion to identify 
the subset of initial values 0(0) and /(0) that lead to asymptotically Lifshitz black 
holes for any value of z > 1. We already learned from (2.19) that the numerical 
integration towards s > 0, i.e. towards the asymptotic region, will terminate at a 
finite value s = s^. We also know that e 2p ^ v ' u " > — > r 2z , as r — > oo for asymptotically 
Lifshitz geometry. For a given 0(0) we tune /(0) until the combination p(s) + z<p{s) 
goes to a finite value in the limit s — > Soo. Furthermore, through the transformation 
rule (2.30) for the conformal factor we can read off the Hawking temperature of the 



4 By the energy argument of [15], a second mode of the linearized system should also be set to 
zero in the z < 2 case, leading to at best a discrete spectrum of Lifshitz black holes at z < 2. 
This condition is too strong, indicating a problem with the background subtraction used in [15]. 
An alternative definition of the energy of asymptotically Lifshitz solutions, based on boundary 
counterterms, recently appeared in [21] for which only a single mode needs to be fine-tuned away 
in order to have finite energy at z < 2. 
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black hole, 



T, 



H 



K 

27 



: lim e-pW-^W, (2.39) 

S-»Soo 



once we have fine-tuned the initial data to give a finite limit. This method works the 
same way for both z > 2 and z < 2. Although the unwanted mode is slowly decaying 
as a function of r when r — > oo in the z < 2 case, it has a growing amplitude as a 
function of s due to the rapid growth of Jf- as s — )■ Sqo. 

A numerical solution for both the ex- 
terior and interior geometry can now be 
obtained by separately integrating towards 
positive and negative s, using the same 
fine-tuned initial data, and then patching 
the two solutions together across s = 0. 
Figure 2 shows a numerical black hole so- 
lution at z = 4 obtained in this manner. 
The curvature singularity is at s — > — oo 
in these coordinates and the numerical 
evaluation will break down before it is 
reached. 

2.3 Global geometry of generic Lif- 
shitz black holes 




Figure 2: Numerical solution for a z = 4 
Lifshitz black hole with p solid, cf> dashed, 
and / dot-dashed. The horizon is at s = 
and r — > oo corresponds to s — > 1. 



We now give an argument that the conformal diagram in Figure 1 applies to a generic 
Lifshitz black hole solution that has a single non-degerate horizon and is obtained 
from the action (1.4), for some value of z > 1. To do that, one examines the 
asymptotic behavior of solutions to the system of equations (2.32) - (2.34) near the 
singularity, i.e. as s — > — oo. The structure of the equations restricts the asymptotic 
behavior to one of two types. The first type has 

-2a 



I: 



:-s) a e*' + ..., 



;_ a )-°- VV° + ..., f = (-s)- m f + ..., (2.40) 



where a, <po, po, and / are non-universal constants. The assumption that r — > 
as s — > — oo amounts to the requirement a > but otherwise these constants are 
unrestricted a priori. 5 

With a > the last term in (2.32), the last two terms in (2.33), and the next to 
last term in (2.34) are sub-leading for this type of solution. The remaining terms in 
(2.32) cancel automatically at leading order, while (2.33) and (2.34) are satisfied at 
leading order if 

f e^ = Z -. (2.41) 



k 



e 2p +2<t>o _ _ Q 2 



5 We note that our argument assumes that r — > as s — > — oo and does not preclude the existence 
of solutions with a smooth inner horzion at which r — > r_ as s — >■ — oo for some finite r_ > 0. On 
the other hand, none of the known exact solutions at z > 1 have such an inner horizon. 
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It follows that this type of asymptotic behavior can only be realized in a k = — 1 
geometry. It is easily checked that the z = 2 topological black hole solution is of this 
type with a — \. 

The other possible behavior near the singularity is 

II: = {-sf e*> + . . . , ef = (-s)- 2/3 " 1/2 e p0 + . . . , / = / + . . . , 

f = {s)-^- 1 h + ..., (2.42) 

with f3 > 0. In this type of solution the last term in (2.32), the fourth and fifth terms 
in (2.33), and the fourth term in (2.34) are subleading as s — > — oo. The remainder 
of (2.32) is automatically satisfied at leading order, while (2.33) and (2.34) require, 

fl e 2p 0+ 4<p = 2> f 2 = f 2 e 2 P0 (2>43) 

This time around there is no restriction on the value of k. The exact z = 4 black 
hole solutions fall into this category, with j3 — § in the k — 1 case and j3 — | in the 
k = — 1 case. 

The Ricci scalar diverges in the s — > — oo limit. For type I solutions the leading 
behavior is R = + ■■■■, while for type II solutions we find R = ^ + . . .. The 
curvature singularity at r — > is always null. This follows immediately from the fact 
that 

e 2p ^0, as s ->■ -oo, (2.44) 

for both allowed types of asymptotic behavior. As a result, the conformal diagram 
in Figure 1 describes generic asymptotically Lifshitz black holes with a single non- 
degenerate horizon that are obtained from the action (1.4). 

In the notation used in [13], appropriately continued to inside the horizon, a 
type I solution has / oc r 1 ~ z , g oc r 2 , h oc 1 and j oc r _1 , for small r. The 
corresponding limiting behavior of a solution of type II is / oc r 2 ~ z , g oc r 2 , h oc r~ 2 
and j oc r~ 2 . 



3. Charged Lifshitz black holes 

We wish to extend the notion of a holographic superconductor to systems with Lif- 
shitz scaling. In order to do this, we need to introduce some additional ingredients to 
the gravity model that we have been considering. As it stands, the model has only a 
single characteristic length scale and does not support a superconducting transition, 
or any other phase transition for that matter. If we want to make phase transitions 
possible, we need to introduce a new scale into the problem. This can be seen explic- 
itly as follows. The gravitational dual of a holographic superconductor is a charged 
plane-symmetric black hole with hair [4] , so for this application we work with k — 
geometries of the form, 

ds 2 = L 2 [-e 2 ^dvdu + e~ 2 ^ (dx 2 + dy 2 )] . (3.1) 
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The plane-symmetric black hole geometry has a scaling symmetry, which is absent 
for spherical or hyperbolic horizons. The metric is invariant under a global rescaling 
of the planar coordinates x, y — > e a x, e a y accompanied by a uniform shift <fi — > <fi + a 
while keeping p, v, and u unchanged. On the other hand, such a shift changes the 
Hawking temperature (2.39). Since different non-zero Hawking temperatures can be 
mapped into each other by scale transformations they are all physically equivalent 
in this system and there cannot be any phase transitions. Another way to state this 
is that in a scale invariant theory at finite temperature there is no other length scale 
apart from that defined by the temperature itself and therefore different temperatures 
must be physically equivalent. 

There are different ways to add another length scale to the Lifshitz system. We 
choose to do it by letting our black holes carry an electric charge, which couples to 
a Maxwell field J- {2)- This corresponds to introducing a non-zero charge density into 
the dual field theoretic system [4] and adds a term to the action, 

S ? = -\f *^(2)A7( 2 ), (3.2) 

which has the same form as the kinetic term of F( 2 ) and contributes in same way to 
the field equations. Introducing a new gauge field may, at first sight, appear to be an 
unnecessary complication given that the model already has a two-form field strength, 
F( 2 ), and a three-form field strength, #(3), which couples to the two- form field and 
could be viewed as a form of charged matter. While such an interpretation could in 
principle be pursued, we find it more useful to consider Ft 2 ) and as auxiliary 
fields, that only couple to gravity, and whose only role is to give rise to gravitational 
backgrounds with Lifshitz scaling. We are then free to separately include the phys- 
ical electromagnetic field J r ( 2 ) and introduce charged matter that couples to it but 
not to the Lifshitz form-fields. The manner in which we construct our holographic 
superconductors with Lifshitz scaling is analogous to the corresponding construction 
for the z = 1 case in [4]. In fact, in the limit z — > 1, our system reduces to the one 
considered by those authors and we have used this to check our numerical algorithm 
against known results for z — 1. 
Writing 

J-"( 2 ) = Lp vu (v,u) dv A du, (3.3) 
the Maxwell equations for J 7 ^) reduce to 

d v {e- 2 ^ Pvu ) =0 = d u (e"^ 2 V) . (3.4) 

The general solution can be written 

Pvu = Qe 2 ^, (3.5) 

and has a simple interpretation as the Coulomb field of a point charge. Now consider 
a black hole with charge Q in the s variable. Since the new gauge field does not 
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couple directly to the original Lifshitz gauge fields, the field equation for / (2.32) 
remains unchanged while equations (2.33) and (2.34) pick up terms involving the 
black hole charge, 6 

= S(f) " + (/,'- 2s 0' 2 + ^ e 2p+2<i> + (^+ 4 ) e 2p _ l (/ ~2 + g2 )e 2 P+ 4^ (3 6) 

= sp " + p '- s 0' 2 + k -e^ - Up + Q 2 )e 2p+i4> - fe^f' 2 . (3.7) 

4 2 2z 

The field equations can be numerically integrated as before and we will present some 
numerical results below, but before that we present some exact solutions. In the first 
example we recover the well known AdS-Reissner-Nordstrom solution as a special 
case with z — 1 and / = 0. This provides a check of the formalism. The remaining 
examples are new and describe one-parameter families of charged Lifshitz black holes 
at z — 4. 

3.1 AdS-RN black holes at z = 1 

The AdS-Reissner-Nordstrom solution at z — 1 describes an electrically charged 
black hole in asymptotically AdS spacetime. In the variables we are using it is given 
by 

= 2p(r») _ r 2 j_ j, _ ±_ I 3 

r \ h r h J r z 

with k = +1,0,-1 for a spherical, flat, or hyperbolic horizon at r = tv The 
relationship between the area and tortoise coordinates is 



e 2p(r,) =r 2 + k _}_( r 3 + krh+ Q^ + Q^^ (g^ 



dr 2,1/,, Q 2 \ Q 2 
dr* r \ r h J r z 



and the Lifshitz gauge field / is everywhere vanishing. It is easily checked that 

e 2p(s) = 1 e 2p(r,) e -0( S ) = (3 1Q ) 

with s = e 2K ( r *- r ? : ) i s a solution of equations (3.6) and (3.7) with / = 0. 
The Hawking temperature, 



T h = — [3r h + --\ , (3.11) 



— 3r h H 3 

4tt V '"h 

goes to zero as the charge approaches the extremal value for a given horizon area, 

T H ^0 as Q 2 ^Q 2 cxt = 3r 4 h + kr 2 h , (3.12) 

as shown in Figure 3. 



6 For completeness we allow for k ^ in the field equations but we are primarily interested in 
the k = case for the application to holographic superconductors. 
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3.2 Exact solutions for charged black holes at z = 4 

Numerical work is required in order to explore the full parameter range of asymp- 
totically Lifshitz black hole solutions at z > 1 but exact solutions are useful, even 
if they only apply in special cases. We have found a family of exact charged black 
hole solutions for z — 4, which generalize the isolated z = 4 solutions discussed in 
Section 2.1. In the notation of equations (2.2) - (2.3) the metric is given by 



/(r) = ^ = V 1+ i^-I5o^-2^' (3 - 13) 

with k = +1, 0, —1. The metric reduces to the z = 4 solutions of Section 2.1 when 
Q — and k — ±1 while the Q = 0, /c = case reduces to the z = 4 Lifshitz fixed 
point geometry. 

The tortoise coordinate is given by equation (2.13) with 



* = Ws> + t-4 (3 ' 14) 



^+«! + * (3.15) 
100 2 20 V ; 

The change to U, V coordinates in (2.14) leads to a metric of the form (2.15), with 
the new values of b\ and which is nonsingular at the horizon. The Lifshitz gauge 
field can be obtained by inserting the p(s) and 0(s) read off from (3.13) into (3.6) 
and solving for f(s). One finds that / = (r 2 + ^) is independent of the black 
hole charge Q. 

The Hawking temperature is Th — ^ with 



The k = metric has 61 = b 2 and takes a particularly simple form in U, V coordinates, 
ds 2 = L 2 (^-^dUdV + r 2 (d6 2 + 9 2 d<p 2 ) S j . (3.17) 

In this case the horizon is at r = ^ and the Hawking temperature is T H — both 
of which go to zero in the Q — > limit. 

Returning to the more general exact solution in (3.13), the Ricci scalar is given 

by 

1 f3Q 2 9k 2 3k 



R = TA^ + w^-^- hA )' (318) 

and is singular at r — > 0. The curvature singularity is null and the conformal diagram 
remains the same as in Figure 1. In particular, unlike ordinary Reissner- Nordstrom 
black holes, these charged Lifshitz black holes do not have an inner horizon away from 
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the singularity at r — > 0. The general analysis from Section 2.3 of the asymptotic 
behavior near the singularity can be extended to the case of charged Lifshitz black 
holes. They exhibit type II behavior with the replacement /q — » /q + Q 2 in equation 
(2.42). 

3.3 Numerical charged black hole solutions 

The field equations with the additional Maxwell field included were presented earlier 
in equations (2.32), (3.6), and (3.7). The initial values giving a smooth horizon at 
s = are modified as follows by a non-vanishing black hole charge, 

0'( O ) = i (-ke 2 ^ - (z2+ 2 Z+4) + (/(O) 2 + Q 2 )e 40(O) ) e 2 ' (0) , (3-19) 
p'(0) = \ (y- k -e 2 ^ + (/(O) 2 + Q 2 )e^ e**°\ (3.20) 
/*(0) = f/(0)e^°). (3.21) 

For any given value of z > 1, there is a two-parameter family of asymptotically 
Lifshitz black hole solutions. If we take the charge Q and 0(0) as independent 
parameters then /(0) needs to be fine-tuned in order to obtain the correct asymptotic 
behavior as s — > Sqo. Changing the value of p(0) amounts to a global rescaling of the 
s coordinate and does not change the geometry. 

The numerical integration proceeds 
in the same way as before and plots of 
numerical solutions are qualitatively sim- 
ilar to Figure 2. The existence of a new 
length scale in the system can be seen in 
the Hawking temperature obtained from 
(2.39) for black holes with different Q 
keeping the area coordinate at the hori- 
zon fixed. Figure 3 shows the black hole 
temperature function of Q for sev- 
eral values of z. All the black holes have 
0(0) = 0. In order to facilitate com- 
parison, the temperature of all the black 
holes for a given z is normalized to the 
Hawking temperature of a Q = black 
hole with the same value of z, and simi- 
larly the charge is normalized to the ex- 
tremal charge of black holes with the given f rom numer i ca i solutions 
value of z. 



T H /T 




0.2 0.4 0.6 0.8 1.0 



Figure 3: The Hawking temperature of 
charged black holes with r^=l, k=+l for 
z = 4, 2, |, 1 from bottom to top. The tem- 
perature is normalized to the temperature 
of the corresponding uncharged black hole 
while the charge is normalized to the cor- 
responding extremal charge. The curve for 
2=1 is plotted from the analytic expression 
(3.11) while the z>l curves are obtained 
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4. Coupling to charged matter 



The final ingredient is a scalar field if), which is charged under the new gauge field 
J r (2) but neutral under the original Lifshitz fields F( 2 ) and H^. This allows our 
black holes to grow scalar hair, and, under certain conditions, the black hole hair 
corresponds to the condensation of a charged operator at low temperatures in the 
dual field theory. The overall picture is quite similar to the one obtained for charged 
AdS black holes in [4] and extends the notion of a holographic superconductor to 
systems that exhibit Lifshitz scaling with z > 1. 

4.1 Charged scalar field 

The scalar field action is given by 

S* = -\J d A x^ (<T(^* + iqA^*)(d„il> - iqA^) + m 2 ^) , (4.1) 

where A^ are the components of a one-form potential for the two-form field strength 
J r ( 2 ) and q is the charge of ip. There is an analog in this system of the Breitenlohner- 
Freedman bound [5] on the allowed mass of the scalar field, 

L 2 m2> _(i+^)! (42) 

4 

The requirement that the corresponding Euclidean action be finite places restrictions 
on the asymptotic behavior of the scalar field as r — > oo. This was worked out in 
detail for the z = 2 case in [12], 7 and, since the analysis carries over to general 
z > 1 in a straightforward way, we merely summarize some results without going 
into details. 

The scalar field equation obtained from (4.1) has two independent solutions 
vp(x^) = c + ip + (x^) + c_'0-(a;' 1 ), that have the asymptotic form 

^ ± ( x »)^ r - A ±4>±(T,e,<p) + ..., (4.3) 

at large r, where 

A± = ZJ r ± \l{^r) 2+m2L2 - (44) 

If the mass squared satisfies 



7 The discussion in [12] was in turn based on earlier work on symmetry breaking in the context 
of the AdS/CFT correspondence [22]. 
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then only ib + falls off sufficiently rapidly as r — > oo and the scalar field has the 
boundary condition 

^af) r -A + ^ +(T) 9j v?) + { 1)^ . (4.6) 

The scalar field is then dual to an operator of dimension A + > | + 2. 
If, on the other hand, the mass squared is in the range 

1 _r£^v>^>-r^v. (4.7) 



then both ■0+ and ■?/>_ have sufficiently rapid falloff and there is a choice of two 
different quantizations for the scalar field. In one case the scalar field is asymptotic 
to tb + and dual to an operator of dimension A + with | + 1 < A + < | + 2. In the 
other case the scalar field is asymptotic to ib- an d dual to an operator of dimension 



A_ with | < A_ < | + 1. 



In our numerical calculations we set the scalar mass squared to 



■2 2 1 /^ + 2 



LW = i -^^j 2 , (4.8) 

which is inside the range where there is a choice of two boundary theories. This value 
leads to convenient values for the operator dimensions, A± = ± \. Non-linear 
descendants of the leading scalar field modes are suppressed by O(^) at r — > oo and 
this choice of mass squared ensures that the first descendant of ■0- falls off faster 
than ib + . 

4.2 Lifshitz black holes with scalar hair 

We will look for static spherically symmetric black hole solutions using a metric of 
the form (2.18). We make the gauge choice, 

A = L(a v dv + a u du) , (4.9) 

with a v = a u = a. In this gauge the Maxwell equations for static configurations 
imply the equation, 

„ dib dib* 

which implies, in turn, that the phase of ib is constant and we can take the scalar 
field to be real valued. The remaining non-trivial Maxwell equation is 

d ( -in-u, da \ 



dr„ 



^ e -2 P -24>d^J + q 2 L 2^ e ~2* a = 0, (4.H) 



which agrees with (3.4) if q = 0, i.e. when ib is a neutral scalar field that does not 
act as a source for the Maxwell field. 
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The next step is to write equations for static configurations using the s variable, 
in which the metric is non-degenerate at the horizon. Since a is a component of a 
one-form potential, it transforms under a change of coordinates, 

ds 

a(r*) = —a(s) = 2K,sa(s). (4-12) 

(XT* 

It follows that the potential vanishes at the horizon in tortoise coordinates as long as 
it remains a smooth function there in the s variable. The Maxwell equation becomes 

and for a smooth horizon at s = the initial data must satisfy 

It 2 

-o'(0) + a(0)(p'(0) + 0'(O)) = i_V;(0) 2 a(0)e 2 ^. (4.14) 

8 

The value of a(0) is a free parameter that determines the black hole charge. To see 
this, we note that for q = we have a + sa' = Q-e 2 ^ +2p , with Q the black hole charge, 
and in this case 

a ( ) = |e 2 ^ +2 ^. (4.15) 

If there is a no n- vanishing scalar field with q ^ 0, then the solution of the Maxwell 
equation will no longer be a simple Coulomb field but we are nevertheless free to 
write a(0) the same way as in (4.15) and the constant Q can still be interpreted as 
the total charge inside the black hole. 

The scalar field equation in the s variable is 

2 T 2 

= s^" + ip' - 2s0V + 4sg 2 L 2 a 2 ^ - ^^e 2 >, (4. 16) 

with the following condition on initial data at a smooth horizon, 

^'(0) = !^L e 2p(o)^ (0) . (417) 

The equation for / remains unchanged but the remaining two field equations get 
contributions from the scalar action, 

= S( p" + 4'- 2S 0' 2 + -e 2 ^ + (^+ 4 ) e 2p _ I / ~2 e 2p+40 

4 8 4 

-4(a+ sa') 2 e- 2 ' - ^^e 2 V, (4.18) 

= sp " + p >- S( p' 2 + -e 2 P^ - ±p e 2p+V> _ ± e 4<t> f 2 

4 2 2^ 

-8(a + s a') 2 e~ 2p + -^' 2 - sg 2 L W, (4.19) 
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and the conditions on the initial data become, 



0'( O ) = 1 ^_^(0)_ (^+^+ 4 ) +( J (0) 2 + g 2 )e 4^(0) + ^^ (0)2 ^ e 2p(0) j(42()) 

p'(0) = \ (-|e 2 ^ + (/(O) 2 + Q 2 )e^A e 2 ^. (4.21) 

For given values of z, k, q, and m 2 there is now a three-parameter family of smooth 
initial values for the dynamical fields at the horizon given by Q, ip(0), and /(0). 
As before, the initial data must be fine-tuned to obtain an asymptotically Lifshitz 
geometry and we take the black hole charge Q and the value of the scalar field at 
the horizon t/j(0) as independent parameters. 



\cJ/T c 




Figure 4: Holographic superconductivity at z = 3/2. The graph shows the condensate in 
the O- theory as a function of temperature normalized to the critical temperature. 

4.3 Superconducting phase 

In order to study black hole geometries that are dual to a holographic superconductor 
we set k = 0. For a given value of z we set the scalar mass squared to the value in 
(4.8) and select a value for the scalar field charge q. We then find numerical black 
hole solutions with scalar hair for a range of Q and ip(0) and study the asymptotic 
large r behavior of the scalar hair in each case. The signal of a superconducting 
condensate in one or the other boundary theory is to have either 

c+ = and (£L) = c_ + 0, (4.22) 

or 

c_ = and (0+) = c+ ^ 0. (4.23) 

The curve of vanishing c + (c_) in the Q vs. tp(0) plane can be tracked, tabulating the 
value of c_ (c+) as a function of the black hole temperature along the way. Figure 4 
shows the result of this procedure for hairy Lifshitz black holes at z = 3/2 for scalar 
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field charge q = 1. There is clear evidence of condensation in the corresponding 
boundary theory. Analogous results are obtained for other values of q and more 
general z > 1 but we defer a more detailed study of the parameter space to [23]. 

5. Conclusions 

In this paper we have given a holographic description of superconductors with Lif- 
shitz scaling. The key ingredient is a family of black holes with scalar hair in an 
asymptotically Lifshitz space time. These black holes are interesting in their own 
right, and we have studied their geometry in detail uncovering several intriguing 
properties, including the presence of null singularities. Our work is an extension of 
previous work on Lifshitz black holes in [13, 14, 15]. 

As in the case of the more conventional holographic superconductors at z — 1, see 
[1, 2], our black holes need to be equipped with additional charges in order for a phase 
transition to a superconducting phase to occur. We provide such a construction, and 
show how a charged scalar in this background condenses at a critical temperature, 
indicating a superconducting phase. 

Our results demonstrate that phase transitions leading to superconductivity at 
low temperature, can be given a holographic description also for Lifshitz points. Such 
superconductors are known from solid state physics, see, e.g. [24], and we hope that 
our results will be useful for further studies of these systems. Some of their properties 
are different from usual superconductors, and it would be interesting to see whether 
this has a holographic counterpart. 
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